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A  Simple  Method  for  Predicting  RF  Attenuation 
Through  a  Rocket  Exhaust  Plume 

L.  Glen  McMillion 
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=  plume  drag  coefficient 

=  exhaust  gas  mass  fraction  of  total  gas  in  the  mixing  region 
=  plume  inviscid  core  length 
=  mass 

=  Mach  number 
=  gas  molecular  weight 

=  gas  phase  chemical  species  number  density,  number/cm3 
=  pressure 
=  plume  radius 
=  gas  constant 
=  Schmidt  number 
=  temperature 
=  velocity 

=  axial  (longitudinal)  position 
=  thickness  of  plasma  slab  stream 
=  gas  specific  heat  ratio  (Cp/Cv) 

=  electron-neutral  body  collision  frequency 
=  gas  density 


=  centerline 
=  nozzle  exit 
=  electron 

=  inviscid  core  boundary  (when  used  for  plume  radius,  ‘r’) 

=  core  value  if  x  <  L;  center  line  value  if  x  >  L  (when  used  for  ‘u’  and  ‘f ) 
=  index  counter  on  summations  (when  used  with  2) 

=  start  of  pressure-balanced  jet  (mixing/afterbuming  region) 

=  sound 
=  turbulent 
=  free  stream 
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1.0  INTRODUCTION 

A  theoretical  description  of  the  governing  mechanisms  of  attenuation  of  radio  frequency 
electromagnetic  waves  through  rocket  motor  exhaust  plumes  is  developed  and  presented  as  a 
simplified  mathematical  model  for  prediction  of  attenuation.  This  presentation  is  suitable  as  the 
basis  for  writing  a  computer  code  to  perform  RF  signal  attenuation  predictions. 

The  existing  theory  for  predicting  radar  attenuation  by  absorption  applies  only  to  weakly- 
ionized  plasma  with  uniform  properties;  the  plasma  resulting  in  the  exhaust  of  a  typical  tactical 
rocket  motor  can  safely  be  assumed  to  be  weakly  ionized,  but  cannot  be  assumed  to  be  uniform 
throughout.  Therefore  we  must  calculate  the  properties  of  several  adjacent  streams,  which  are 
referred  to  as  slabs,  within  the  exhaust  plume,  each  having  approximately  uniform  properties,  and 
transversed  by  the  radar  beam  in  a  series  sequence.  The  total  attenuation  is  the  summation  of 
transmission  loss  through  the  adjacent  slabs. 

Therefore  predicting  radio  frequency  transmission  loss  through  a  rocket  exhaust  plume 
requires  that  the  exhaust  plume  model  provide  a  description  of  the  spatial  variation  of  temperature 
and  chemical  species  concentrations.  Models  of  varying  degrees  of  complexity  exist.  The  model 
presented  here  represents  the  least  complex  approach  that  has  been  found  to  provide  reasonably 
accurate  results.  This  model  consists  of  the  following  steps: 

1 .  Calculate  nozzle  exhaust  conditions. 

2.  Calculate  correction  for  non-optimum  nozzle  expansion. 

3.  Calculate  plume  inviscid  core  flow  field. 

4.  Calculate  plume  mixing/afterbuming  region  flow  field. 

5.  Calculate  plume  electron  and  neutral  body  density  field. 

6.  Calculate  electron-neutral  body  collision  frequency. 

7.  Calculate  radio  frequency  transmission  loss. 
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2.0  TECHNICAL  DISCUSSION 

2.1  Nozzle  Exhaust  Conditions 

Nozzle  exhaust  conditions  are  calculated  using  the  Naval  Weapons  Center  Propellant 
Evaluation  Program.  There  are  two  versions:  MicroPEP  and  PEP93.  PEP93  is  a  newer  version 
and  is  more  user-friendly  than  MicroPEP.  However,  the  input  format  for  MicroPEP  may  make  it 
easier  to  incorporate  into  a  larger  exhaust  plume  RF  interaction  code.  Copies  of  both  programs 
are  included  with  this  report. 

Input  data  required  to  run  PEP  to  calculate  nozzle  exhaust  conditions  are: 

1 .  Propellant  formulation.  The  propellant  formulation  MUST  include  concentrations  of  the 
alkali  metals  sodium  and  potassium.  The  attenuation  of  RF  signals  by  rocket  exhausts  is  due  to 
free  electrons  in  the  plume,  which  are  formed  primarily  by  ionization  of  sodium  and  potassium. 
Sodium  and  potassium  are  not  usually  present  as  intentional  ingredients  but  are  present  as 
contaminants  in  other  ingredients,  particularly  ammonium  perchlorate.  Each  ingredient  in  the 
formulation  must  be  assayed  for  sodium  and  potassium  and  the  total  amount  of  each  calculated 
for  the  propellant.  Any  sodium  and  potassium  present  as  contaminants  will  be  present  as 
positively  charged  ions,  not  as  the  metal.  Therefore  an  ionic  species  such  as  sodium  chloride  or 
potassium  chloride  should  be  chosen  from  the  propellant  ingredient  list.  The  weight  of  each  salt 
to  use  must  be  adjusted  to  account  for  the  weight  of  the  anion  (chloride,  in  this  example)  while 
providing  the  total  weight  of  alkali  metal  present  as  calculated  from  ingredient  assays.  Note  that 
the  amount  of  alkali  metal  salts  input  as  propellant  ingredients  must  include  the  contribution  from 
nozzle  throat  and  exit  cone  ablation,  which  is  discussed  next. 

2.  Alkali  metal  mass  flow  rate  during  motor  operation  due  to  ablation  of  nozzle  throat  and 
exit  cone.  Nozzle  throats  and  exit  cones  are  usually  made  of  materials  that  are  ablative  and  that 
contain  sodium  and  potassium  as  trace  contaminants.  The  ablation  rates  must  be  known  so  that 
the  mass  flow  rate  of  alkali  metals  introduced  into  the  plume  by  nozzle  ablation  can  be  calculated. 
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Ablation  rates  are  usually  given  as  millimeters  per  second  or  a  similar  length  per  unit  time  value. 
To  calculate  the  mass  flow  rate  of  ablative  material  the  area  of  ablative  material  exposed  to  the 
combustion  gases  must  be  calculated  for  the  nozzle  geometry.  Alkali  metal  mass  flow  rate  is  then 
calculated  by  multiplying  the  ablative  material  mass  flow  rate  by  the  mass  fraction  of  alkali  metal 
in  the  ablative  material.  The  total  amount  of  alkali  metals  introduced  into  the  plume  during  the 
motor  bum  is  found  by  multiplying  the  mass  flow  rate  by  the  bum  time.  This  number  can  be 
treated  as  an  equivalent  amount  of  alkali  metals  in  the  propellant.  This  is  how  alkali  metals  due  to 
nozzle  ablation  should  be  included  in  the  calculation. 

3.  Rocket  chamber  operating  pressure. 

4.  Nozzle  expansion  ratio. 

5.  Nozzle  exit  pressure  (ambient  pressure). 

When  running  PEP  the  following  options  must  be  selected:  the  program  must  be 
instructed  to  include  ionic  species  in  the  calculations,  the  lower  limit  of  species  concentrations 
must  be  increased  by  at  least  four  orders  of  magnitude  (the  default  limit  is  10"6),  the  density 
exponent  should  be  unity,  and  all  combustion  species  used  in  the  calculation  should  be  listed. 

Chemical  reactions  are  considered  to  be  frozen  through  the  nozzle  expansion  rather  than 
assuming  that  chemical  equilibrium  is  maintained  during  the  expansion  process.  It  has  been  found 
that  equilibrium  chemistry  calculations  applied  to  nozzle  expansion  under-predict  trace  species 
concentrations,  especially  free  electron  concentration,  that  are  critical  to  radar  attenuation 
predictions. 

2.2  Nozzle  Non-Optimal  Expansion  Correction 

When  a  rocket  exhaust  nozzle  is  not  optimally  expanded  resulting  in  the  nozzle  exit 
pressure  being  not  equal  to  the  ambient  pressure,  a  correction  is  needed  to  locate  the  axial 
distance,  Xj,  at  which  pressure  equilibration  is  attained  and  to  determine  the  maximum  or  minimum 
inviscid  plume  radius,  which  occurs  at  the  axial  location  Xj. 
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For  low  altitude  plumes,  if  the  jet  properties  (re,  Te,  Pe)  are  known  at  the  rocket  nozzle  exit 
(calculated  here  from  the  NWC  MicroPEP  code),  these  may  be  extrapolated  to  balanced  pressure 
conditions  (Pj  =  P„)  assuming  conservation  of  mass,  momentum,  and  energy.  In  this 
extrapolation,  no  mass  entrainment  is  allowed,  so  that  mixing  with  ambient  or  free  stream  air 
starts  only  after  pressure  equilibration  is  reached. 

It  is  reasonable  to  assume  that  for  the  purposes  of  mixing  model  definition,  the  balanced 
pressure  is  reached  at  about  the  middle  of  the  first  plume  wavelength.  This  may  be  obtained  from 
(Jarvenin  and  Hill,  1970): 

Xj  =  0.69  re  Me  (ye  PJPJ  0  5  (1) 

At  this  axial  position,  plume  velocity,  temperature,  and  radius  may  be  given  by  (Sukanek,  1976): 

Vj  =  Ve  0  +  0  (2) 

T;  -  [1  -  (Y,  -  1)  5  Ml  -  0.5  (Y,  -  1)  e  Ml]  Te  P) 


\p 

T  v 

e 

j  e 

T  v 

where: 


-B  +  \jB2  +  4 AC 
2A 


(5) 


A  =  Ye  M2e  +  p 


(6) 


B  =  (y,  -  1)  P  M2e  -  a  +  ye  M2e 


(7) 
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C  =  a  +  P 


(8) 


(9) 


P  =  0.5  Y-  CD  Ml 


(10) 


For  low  altitude  plumes  RF  attenuation  is  dominated  by  the  near-field  plume,  that  is,  the  inviscid 
and  mixing/afterbuming  regions  of  the  plume.  For  this  case  it  is  reasonable  to  neglect  plume  drag 
resulting  in  CD  =  0  and,  consequently,  P  =  0.  This  assumption  simplifies  the  equations  to: 


a  =  Ye  Ml 

B  =  Ye  Ml  -  « 

C  =  a 

1  \P~ 
a  =  1  -  — 

Pe\ 

P  =0 


(11) 

(12) 

(13) 

(14) 

(15) 


Equation  (1)  is  used  to  find  the  axial  position,  Xj,  at  which  the  pressure-balanced  jet 
begins.  Equation  (4)  is  used  to  find  the  jet  radius,  rJ5  at  Xj.  The  jet  radius  rj  is  used  as  the  starting 
radius  for  subsequent  calculations  of  jet  properties.  Equations  (2),  (3),  (5),  and  (11)  to  (15)  are 
used  to  solve  equation  (4)  for  r^ 
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2.3  Plume  Inviscid  Core  Flow  Field 

The  procedure  given  in  Section  2  above  is  a  simplified  means  of  allowing  for  under-  or 
over-expansion  of  the  rocket  nozzle.  Its  main  effect  is  to  change  the  scale  of  the  jet;  the  part  of 
the  jet  downstream  of  Xj  behaves  as  if  it  were  produced  by  a  parallel-flow  nozzle  of  radius  r^.  The 
inviscid  core  and  mixing  layer  scales  are  accordingly  normalized  by  reference  to  the  radius  rj 
rather  than  re. 

The  inviscid  core  length,  L,  is  assumed  to  be  a  function  of  the  Mach  number,  Mp  on  the  jet 
center  line  at  the  axial  position  Xj,  and  is  given  by: 

-  =  2.1  M2  (16) 

rJ 

where  ij  is  calculated  from  equation  (4). 

The  Mach  number,  Mj,  is  the  ratio  of  the  local  gas  velocity  to  the  local  speed  of  sound; 

M  =  i  (17) 

1  v 

The  speed  of  sound,  vs,  at  Xj  must  be  calculated  from  the  gas  properties  at  xf 


VPJ0'5 
P V  . 


(18) 


The  gas  density,  pj?  may  be  calculated  using  the  ideal  gas  law,  and  the  gas  temperature,  pressure, 
and  molecular  weight  (MW)  at  Xj : 


Pj 


=  MW 


R  T 


(19) 
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Substituting  equation  (19)  into  equation  (18)  yields: 
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Yj  P„  R  z;.]0-5 
MW  Pj 


(20) 


However,  since  at  Xj  the  jet  pressure  is  equal  to  the  ambient  pressure  (Pj  -  PJ,  equation  (20) 
simplifies  to: 


v  = 


r /* 


MW 


0.5 


(21) 


Now  the  Mach  number  at  jq  can  be  calculated  from  equations  (2),  (3),  (17),  and  (21),  given  y,  and 
MW  from  MicroPEP  calculations.  Gas  molecular  weight  and  y  can  be  considered  constant 
throughout  the  inviscid  core  without  significant  error,  therefore  MW  and  y  from  the  MicroPEP 
exit  plane  calculations  will  be  used.  Once  Mj  is  known  the  core  length  L  can  be  calculated  from 
equation  (16). 

The  inviscid  core  is  assumed  to  be  a  cone  of  base  radius  r}  and  height  L,  so  that  the 
generator  is  given  by: 


(22) 


2.4  Plume  Mixing/Afterburning  Region  Velocity  Flow  Field 

The  gas  velocity  is  assumed  to  remain  constant  throughout  the  core,  and  is  allowed  decay 
on  the  axis  beyond  the  core  inversely  with  x: 

L  (23) 

JC 


u 


CL 


8 


V  • 


if 


The  radial  variation  of  the  longitudinal  component  of  velocity  in  the  mixing  region  is 
expressed  as  an  exponential  function  of  position: 


u 

—  =  exp 


(24) 


where: 


Ui  =  iw  if  x  <  L; 

U;  =  uCL  if  x  >  L ; 

r(  =0  if  x  >  L ; 

u. 

rm  =  the  radial  coordinate  at  which  u  =  — 


It  can  he  shown  that: 


where: 


",  ri 


2 a  2 a2 


(25) 


a2  =  - 1WQ -  (26) 

(  rm  ~  ri  )  2 

Equation  (25)  is  a  quadratic  in  ‘a’  as  given  by  equation  (26)  and  therefore  in  rm .  If  x  <  L  , 
then  u.  =  u e  and  rm  may  be  expressed  in  terms  of  r  ,  which  is  a  fiinction  of  x  as  given  by 
equation  (22).  If  x  >  L  then  r.  =  0  and  =  ue  (L  I  x)  so  that  rm  may  be  expressed 
directly  as  a  function  of  x.  This  defines  the  complete  plume  flow  field. 


2.5  Plume  Electron  Density  and  Neutral  Body  Field 

A  rigorous  determination  of  the  electron  density  field  requires  knowledge  of  the  gas 
temperature,  chemical  composition  and  chemical  reaction  rates.  This  means  that  an  exact  solution 
would  be  undesirably  complex  in  the  context  of  a  “simple”  model.  Therefore  a  compromise 
between  frozen  and  equilibrium  chemistry  with  respect  to  electron  density  distribution  will  be 
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used.  If  the  chemistry  of  the  plume  downstream  of  the  nozzle  exit  plane  is  assumed  to  be  frozen 
the  plume  must  behave  in  all  respects  as  a  simple  jet  of  heated  inert  gas.  Under  this  assumption 
the  mass  fraction  ‘f  of  initial  exhaust  gas  per  unit  mass  of  total  gas  (exhaust  gas  plus  entrained 
free  stream  air)  in  the  mixing  region  is  described  by  observed  empirical  relationships.  The  mass 
fraction  of  exhaust  gas  to  total  gas  in  the  plume  mixing  region  is  defined  as: 

/  =  - —  (27) 

+  m 

e  oo 

The  exhaust  gas  concentration  on  the  plume  axis  has  been  found  to  vary  as  x  2  and  the  radial 
profile  of  exhaust  gas  concentration  is  governed  by  the  similarity  relation: 

/  _ 

/, 

where: 

f,  =  f core  if  X  *  L> 

fi=fCL  ifx>L 

The  first  step  in  determining  the  electron  density  distribution  in  the  model  is  to  assume  that  the 
mixing  region  geometry  is  defined  by  the  velocity  flow  field  presented  in  the  previous  section  and 
that  the  pre-combustion  composition  of  the  gas  mixture  is  described  by  the  variation  of  T  with 
position.  Equation  (28)  is  assumed  also  to  apply  to  combusting  systems  and  the  turbulent 
Schmidt  number  is  assumed  to  be  unity.  These  simplifications  really  mean  that  we  are  assuming 
instantaneous  mixing  of  the  exhaust  gas  with  air  within  the  confines  of  the  velocity  field.  The 
value  of  cf  at  any  point  is  given  by: 


Sc, 


(28) 


/  =  1  inside  the  core 


(29) 


f  =  f  exp 


(30) 
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¥ 


f,=  1  for  x  <  L  ;  f  =  ~  for  x  >  L. 

x2 

Having  in  this  way  specified  a  mixture  composition  at  all  points  within  the  mixing  region,  the 
effects  of  secondary  combustion  and  heat  release  are  simulated  by  permitting  instantaneous 
reaction  to  chemical  equilibrium. 

Plume  mixing  region  chemical  equilibrium  calculations  are  performed  using  the  NWC 
MicroPEP  or  PEP93  code;  by  using  the  appropriate  switches  PEP  will  perform  chemical 
equilibrium  calculations  of  any  type.  For  MicroPEP  this  option  is  selected  by  entering  a  “1"  in  the 
seventh  space  of  the  fourth  line  of  the  input  file  and  replacing  the  chamber  pressure  with  the 
desired  temperature  on  the  Final  Lines  of  the  input  file  (see  the  “readme.txt”  file  included  on  the 
MicroPEP  diskette).  The  input  for  performing  this  calculation  consists  of  the  chemical 
composition  of  the  exit  plane  exhaust  gas  and  the  temperature  of  the  mixed  plume  gas  (before 
allowing  for  chemical  reaction)  as  determined  by  mixing  the  exhaust  gas  with  ffeestream  air  at  the 
ratio,  T ,  calculated  from  equations  (29),  (30),  or  (31).  This  is  a  simple  trial  and  error  calculation 
which  finds  the  temperature  of  the  mixed  exhaust  gas  and  air  by  matching  the  enthalpy  lost  by  the 
exhaust  gas  to  the  air  with  the  enthalpy  gained  by  the  air  from  the  exhaust  gas  during  the  mixing 
process.  The  constant-pressure  heat  capacity  as  a  function  of  temperature,  Cp(T),  must  be  known 
for  air  and  the  exhaust  gas  over  the  temperature  range  of  interest.  The  heat  capacity  of  the 
exhaust  gas  must  be  calculated  from  the  heat  capacities  of  the  individual  component  species  and 
their  mole  or  mass  fractions,  depending  on  the  units  of  the  Cp  data  used.  The  heat  capacity  of  air 
is  tabulated  in  many  sources  in  the  open  literature.  Once  the  temperature  (before  allowing  for 
chemical  reaction)  and  composition  of  the  mixed  plume  gas  are  known,  PEP  is  run  to  calculate 
equilibrium  composition  and  enthalpy  released  by  the  chemical  reactions.  The  equilibrium 
composition  includes  free  electrons  and  the  neutral  bodies  that  will  be  required  to  calculate 
collision  frequency  in  Section  2.6.  The  final  temperature  of  the  mixed,  reacted  gas  is  calculated 
by  dividing  the  enthalpy  of  reaction  by  the  constant-pressure  heat  capacity  of  the  mixed  gas  to 
find  the  temperature  increase  due  to  reaction  and  adding  this  value  to  the  non-reacted  temperature 
of  the  system. 
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For  this  simple  model,  enthalpy  changes  due  to  changes  in  momentum  (mixing  and 
velocity  changes)  are  considered  to  be  negligible  with  respect  to  static  enthalpy  of  the  exhaust  gas 
and  free  stream  air.  If  evaluation  of  the  use  of  the  model  indicates  that  conservation  of 
momentum  is  significant,  the  relationships  presented  by  A.  C.  Victor,  NWC  TP  5319,  Part  1, 
Section  2.2,  page  1 1-13  should  be  used  (Victor,  1975). 

2.6  Electron-Neutral  Body  Collision  Frequency 

The  electron-neutral  body  collision  frequency  at  any  point  in  the  plume  is  calculated  from 
(Cummings,  Williams  and  Wilson,  1970). 

v  =  V  2",  (N,  ?,)  (32) 

where: 


V 


N 


8  kT 


n  m. 


mean  thermal  electron  velocity,  where. 


k  =  Boltzmann  constant  =  1.38066x70  23  J  K  l 
T  =  temperature  as  degrees  Kelvin  ( K) 
me-  =  electron  mass  =  9.10953x70  28  gram 


(33) 


V  =  number  density  of  species  i 

qi  -  electron  collision  cross-section  of  species  i 

In  practice,  only  a  few  of  the  species  present  in  the  exhaust  gases  have  collision  cross- 
sections  large  enough  to  contribute  significantly  to  the  collision  frequency.  These  are  typically  H2, 
H20,  N2,  CO,  and  C02.  The  contribution  to  collision  frequency  by  all  other  species  is  neglected  in 
this  model.  Thus  the  calculation  of  the  equilibrium  composition  and  temperature  in  the  mixing 
region  serves  as  the  basis  for  determining  collision  frequency. 
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Electron-neutral  body  collision  cross-sections  are  dependent  upon  the  electron  velocity, 
ve  - .  The  electron-velocity-dependent  collision  cross-sections  in  the  c.g.s.  system  of  units  for  H2, 
H20,  N2,  CO,  and  C02  are  given  in  Table  1  (Smoot  and  Underwood,  1966): 

Table  1 

Collision  Cross-Sections  of  Electrons  with  Neutral  Species 


Species 

Cross-Section,  qj 

(c.g.s  units') 

h2. 

1.45  23  v  +  8.9”16 

e 

h2o. 

5.9  v? 

e 

n2. 

3.29~23  v  - 

e 

CO. 

2.08~23  v-  +  2  46~16 

co2. 

4.7  8  v".1 

e 

2.7  Radio  Frequency  Signal  Transmission  Loss 

A  knowledge  of  Ne  ,  v,  and  radar  frequency,  co ,  is  sufficient  to  allow  a  direct 
prediction  of  radar  attenuation.  The  existing  theory  for  predicting  radar  attenuation  by  absorption 
applies  only  to  weakly-ionized  plasma  with  uniform  properties;  the  plasma  resulting  in  the  exhaust 
of  a  typical  tactical  rocket  motor  can  safely  be  assumed  to  be  weakly  ionized,  but  cannot  be 
assumed  to  be  uniform  throughout.  Therefore  we  must  calculate  the  properties  of  several 
adjacent  streams,  which  are  referred  to  as  slabs,  within  the  exhaust  plume,  each  having 
approximately  uniform  properties,  and  transversed  by  the  radar  beam  in  a  series  sequence.  The 
total  attenuation  is  the  summation  of  transmission  loss  through  the  adjacent  slabs,  each  of 
thickness  yt  : 


y> v,  (Ne~\ 


or 


db  =  0.46  S”=1 


(34) 


Equation  (34)  has  been  used  to  predict  X-band  and  K-band  radar  transmission  loss  for  two 
solid  rocket  propellants  (Smoot  and  Underwood,  1966). 

3.0  SUMMARY  AND  RECOMMENDATIONS 

Classified  and  unclassified  literature  was  searched  for  articles  and  reports  about  methods 
for  predicting  radio  frequency  signal  transmission  loss  through  rocket  motor  exhaust  plumes.  The 
resulting  literature  was  examined  and  evaluated  to  create  a  simple  yet  effective  prediction 
methodology.  A  great  deal  of  the  literature  deals  with  very  complex  methods  or  rocket  flight 
profiles  that  are  of  no  interest  to  this  project.  No  single  article  or  report  presented  a  method  that 
could  be  applied  to  this  project;  many  articles  had  to  be  drawn  upon  and  scrutinized  to  create  the 
prediction  method  presented  in  this  report.  The  resulting  method  is  based  upon  simplifying 
assumptions  that  are  specific  to  low  altitude,  low  Mach  number  flight  profiles  and  propellants  that 
are  suitable  for  tactical  missiles.  The  equations  describing  rocket  exhaust  plume  velocity  fields, 
chemical  and  physical  properties,  and  the  interaction  of  radio  frequency  signals  with  exhaust 
plume  plasmas  have  been  described  and  organized  in  a  way  that  can  be  used  to  write  a  computer 
code  to  calculate  line-of-sight  RF  signal  attenuation  through  an  exhaust  plume. 

All  of  the  effort  to  date  has  been  spent  on  analyzing  the  available  methodologies  and 
creating  a  suitable  methodology  for  this  project  based  on  appropriate  simplifying  assumptions 
(Phase  I);  the  work  presented  in  this  report  required  greater  effort  than  anticipated  and  a  second 
phase  (Phase  II)  will  be  required  to  complete  the  project  as  described  in  the  Statement  of  Work. 
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This  report  presents  the  results  of  Phase  I:  the  complete  theoretical  description  required  to 
write  a  computer  code  for  predicting  RF  transmission  loss  through  rocket  exhaust  plumes. 
Computer  programming  for  Phase  II  will  use  FORTRAN  or  ANSI  C  programming  language  at 
the  direction  of  the  customer.  The  estimated  effort  to  perform  Phase  II  of  this  task  is: 


a.  Type  of  Contract: 

b.  Number  of  Days: 

c.  Labor  Rate: 

d.  Cost  of  Labor: 

e.  Cost  of  Materials: 

f.  Total  Estimated  Cost: 


Time  and  Materials 
340 

$310.00/day 
$105,400.00 
$  3,000.00 
$108,400.00 
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